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The dynamics of a classical heavy particle moving in a quantum environment is determined by a 
Langevin equation which encapsulates the effect of the environment-induced reaction forces on the 
particle. For an open quantum system these include a Born-Oppenheimer force, a dissipative force 
and a stochastic force due to shot and thermal noise. Recently it was shown that these forces can be 
expressed in terms of the scattering matrix of the system by considering the classical heavy particle 
as a time-dependent scattering center, allowing to demonstrate interesting features of these forces 
when the system is driven out of equilibrium. At the same time, it is well known that small changes 
in a scattering potential can have a profound impact on a fermionic system due to the Anderson 
orthogonality catastrophe. In this work, by calculating the Loschmidt echo, we relate Anderson 
orthogonality effects with the mesoscopic reaction forces for an environment that can be taken out 
of equilibrium. In particular we show how the decay of the Loschmidt echo is characterized by 
fluctuations and dissipation in the system and discuss different quench protocols. 
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I. INTRODUCTION 

Understanding the effect of fluctuations and dissipa¬ 
tion in non-equilibrium settings is essential for an ulti¬ 
mate control of quantum systems. Dissipation is on one 
hand unavoidable in realistic systems, and known to play 
an important role in their dynamics — a paradigm is 
the exponential supression of quantum tunneling out of 
a metastable state as modeled by Caldeira and LeggetlP^ 
— while non-equilibrium can provide new levels of tun- 
ability. This is a topic of renewed interest in view of cur¬ 
rent experiments which explore the possibility of quan¬ 
tum information processing, by embedding a qubit de¬ 
gree of freedom in a mesoscopic system.l^ The coupling 
of the qubit to an environment causes decoherence and 
consequently loss of information, which is closely related 
to the fluctuations and dissipation in the system.l^ 

In this context, the quantum Loschmidt echo, also 
known as fidelity, is a useful quantity that indicates the 
sensitivity of the system to small perturbations.!^ In its 
generalization to many-body systems,!^ the Loschmidt 
echo corresponds to the off-diagonal element (norm- 
squared) of the reduced density matrix for the qubit de¬ 
gree of freedom, and its decay in time characterizes the 
environment-induced decoherence.!^ For a fermionic envi¬ 
ronment, this decay is directly related to the Anderson 
orthogonality catastrophe, which describes the response 
of the fermionic system to a sudden perturbation.E^lln his 
seminal work,!^ Anderson showed that the many-body 
ground state of a fermionic system is, in the thermody¬ 
namic limit, orthogonal to that of the same system in 
which a local scattering potential is introduced. More 
precisely, the overlap of the two states decays as a power- 
law with the system size, with an orthogonality exponent 
characterized by the scattering phase shift produced by 
the scattering potential. The orthogonality catastrophe 
plays an essential role in describing the so-called “impu¬ 


rity problems” in which a local degree of freedom inter¬ 
acts with a fermionic environment 

A class of impurity problems is that of a “heavy parti¬ 
cle” embedded in a quantum environment. The impurity 
in this case is heavy compared to those particles com¬ 
prising the environment, and can be treated as a classical 
degree of freedom with semiclassical dynamics dictated 
by the back-action of the environment. This dynamics 
can be described in terms of a Langevin equation, which 
is a stochastic equation of motion that describes at an 
effective, macroscopic level, the effects of dissipation and 
fluctuations induced by the environment on the heavy 
particle. An interesting question is how the dynamics of 
the heavy particle and orthogonality effects of its environ¬ 
ment are related. For a quantum environment in equilib¬ 
rium, the Anderson orthogonality catastrophe exponent 
was conjectured by Sols and GuinecP^to be proportional 
to the dissipation coefficient a heavy particle experiences 
when moving in a metallic environment. This relation 
was later proved to be valid, in the small-distance limit,!^ 
for a heavy particle moving in a quantum environment at 
zero temperature.!^ For non-equilibrium fermionic sys¬ 
tems this probl em has been studied in the context of 
concrete models.!^ ^ ! ^^ ! ^^ ! ^^ !^ 

Motivated by these findings, in this work we calculate 
the Loschmidt echo for small changes of a scattering po¬ 
tential, in a fermionic open system which is taken out of 
equilibrium by imposing a voltage bias. With the aim 
of exploring the relation between the orthogonality ex¬ 
ponent and the dissipation coefficient in this case, we 
express the decay in time of the Loschmidt echo in terms 
of the coefficients of the corresponding Langevin equa¬ 
tion - in particular in terms of the dissipation and noise 
coefficients. To this effect we make use of the recent de¬ 
veloped formalism that describes the effective forc es in 
the Langevin equation in terms of scattering theory.!23l2ll 
Our results apply generally to systems for which changes 
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in the scattering potential can be treated pertubatively. 

The manuscript is organized as follows. We start in 
section Sec. [IT] by presenting the Langevin equation in 
terms of current-induced forces, and the associated force- 
force noise correlator. In Sec. Ill we perform a perturba¬ 
tive expansion of the Loschmidt echo and show that it can 
be expressed in terms of the noise correlator, and discuss 
different quench protocols. In Sec. |IV| we make use of the 
results of Secs. [TT| and III to show that, in equilibrium and 
for zero temperature, the decay of the Loschmidt echo is 
a powerlaw with an exponent dictated by the dissipation 
coefficient a heavy particle experiences in the fermionic 
environment, in agreement with the known orthogonality 
results. Finite temperatures however render the decay 
exponential. In Sec. |V] we turn to the non-equilibrium 
case for which we calulate the decay of the Loschmidt 
echo within linear response in the applied bias. In this 
case we show that the decay of the Loschmidt echo cannot 
be expressed solely in terms of the dissipation coefficient, 
providing a general expression for the decay in terms of 
the macroscopic Langevin parameters. We then discuss 
different time scales for which the results can be cast in a 
simple form. In Sec. |VI| we apply our results to a simple 
example and check the limits of validity of our approxi¬ 
mations, while we list our main conclusions in Sec. m 
Quite a few calculations in this work are rather lenghty. 
To improve readability, and at the same time to make 
this paper self contained, we have included some details 
of these calculations in the Appendices. 


II. LANGEVIN EQUATION AND NOISE 
CORRELATOR 

In this section we briefly review the elements of the 
Langevin equation that governs the stochastic dynam¬ 
ics of a heavy particle in an open quantum environment. 
Throughout this paper we will consider a fermionic quan¬ 
tum environment that can be taken out of equilibrium 
by a difference of chemical potential in the leads as illus¬ 
trated in Fig.[^ In a concrete example, the heavy particle 
can represent the classical vibrational degrees of freedom 
of a molecule or suspended carbon nanotube connected 
to conducting leads. The heavy particle is represented 
through classical degrees of freedom which are coupled to 
the quantum environment, and disturb it as they evolve 
in time. The back-action of this disturbance onto the 
heavy particle gives rise to reaction forces,!^ also called 
current-induced forces in a quantum transport setup. In 
the adiabatic limit, for which the dynamics of the heavy 
particle is much slower that that of the quantum environ¬ 
ment, this effect is well described semiclassically at the 
level of a Langevin equation obtained by tracing out the 
quantum environment. If we denote the degrees of free¬ 
dom of the heavy particle by X(t), the Langevin equation 
reads (in what follows we ommit the time dependence for 



FIG. 1. In the scattering region electrons couple X via the 
scattering potential Vx. When the scatterer moves by 5X, the 
scattering potential changes accordingly, i.e. Vx+sx. This 
gives rise to the two different Hamiltonians Hi = Ho + Vx 
and Hi = Wo + Vx+<sx. The scattering region is finite and the 
dwell time to gives the timescale the electrons spend within. 


notational simplicity) 

- F^^(X) = p„(x) - ^ r^0{x)Xp + . (1) 

p 

On the left hand side. Pa denotes the canonical momen¬ 
tum of coordinate Xa (a = 1, ■ • ■, N), and we have in¬ 
cluded the possibility of an external classical force F'^^(X) 
(troughout the text we will indicate matrices and vec¬ 
tors in the space spanned by the Xa with bold let¬ 
ters) . The right hand side of Eq. 0 contains the forces 
due to the quantum environment. F(X) is the usual 
Born-Oppenheimer force, while the symmetric and an¬ 
tisymmetric parts of the tensor r(X) represent a dissi¬ 
pative and Lorentz-like (and therefore non-dissipative) 
force, respectively. Fluctuations due to shot and thermal 
noise are taken into account by the stochastic Langevin 
force ^(X). 

When the fermionic system is taken out of equilib¬ 
rium, the current-induced forces present qualita tive d if- 
ferences with respect to the equilibrium situation.l^^^^SMl 
The Born-Oppenheimer force F is non-conservative in 
this case, and therefore provides a way of exchanging en¬ 
ergy between the classical field and the quantum environ¬ 
ment which is non-dissipative. The tensor F constitutes 
the first order correction in an adiabatic expansion to 
the Born-Oppenheimer force. It can be split into sym¬ 
metric and antisymmetric components. The antisym¬ 
metric component is a Lorentz-like term, which can be 
interpreted as an effective magnetic field acting on the 
space spanned by X — this term is not relevant for the 
Loschmidt echo (which involves only symmetric compo¬ 
nents as we will se below) and hence will be not dealt 
further with within this work. We denote the symmetric 
dissipative term of F by 7. It is convenient to express the 
latter as 7 = 7°*! -I- 7'^®'!, where 7“'i represents a pure 
non-equilibrium contribution while 7®'^ is a straightfor¬ 
ward generalization of the equilibrium contribution eval¬ 
uated in an non-equilibrium environment. Explicit ex¬ 
pressions for these quantities are given in App. Ac¬ 
cording to our definition, 7®'^ connects to the equilibrium 
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results, but it also contains non-equilibrium terms for fi¬ 
nite bias. The pure non-equilibrium term can take 
negative values and, moreover, render the full dissipative 
term negative. 

Of particular importance for the following discussion is 
the stochastic component characterized by the force- 
force noise correlator. We will see in the next section 
that the Loschmidt echo can be, perturbatively for small 
displacements SX, written in terms of the noise correlator 

i^a/3(t,0 = {(CaW6(0>},, (2) 

where = {Map -f Mpa) /2 indicates the symmet¬ 

ric component of a generic matrix M. To give a con¬ 
crete expression for this noise correlator we consider a 
generic, albeit non-interacting many-body Hamiltonian 
which depends parametrically on time via the potential 
b’xj = 'Ho + Vx- This potential represents the cou¬ 
pling between the heavy particle and the fermionic en¬ 
vironment. The current-induced force operator is given 
by 


^{t) = -VxHx(t). ( 3 ) 

The Langevin Eq. 0 is obtained by calculating the 
quantum-statistical average {^{t)) within an adiabatic 
expansion to linear order in the velocity X, together with 
the quantum and thermal fluctuations given by 

m = T{t) - {T{t)) . ( 4 ) 

The coefhcients of this expansion are instantaneous: 
the noise is assumed to be delta correlated D(t, <') —>• 
D(X)(5(t — t'), and there is no retardation kernel for 
the dissipative term 7 (X). These are the zero-frequency 
limit, respectively, of the force-force correlator 0 and 
force susceptibility 

xrp{t.t') = -^e{t-t'){[Fa{t). Fp{t')]) ( 5 ) 

where 9{t) is the usual step function. In equilibrium and 
assuming steady state, so the time dependence in the 
relevant quantities is through the time difference (t — t'), 
we have 


i7p{^) = - 


ImXajH 


( 6 ) 


where denotes the Fourier transform of the fric¬ 
tion Kernel t').I^The correlator of the fluctuating 

Langevin force and the friction tensor are related in equi¬ 
librium via the finite frequency fluctuation-dissipation 
theorerrPi^ 


D(a;) = w coth (7) 

(we take the Boltzmann constant ks = 1), where D(a;) is 
the real part of the Fourier transform of the fluctuating 
force correlator in Eq. ([^. In the limit w ^ T Eq. Q 


reduces to the classical identity 

Del = 2T7«^ (8) 

where 7 ®'^ is evaluated at zero frequency, while for w ^ T 


D» = |w| 7 '^». ( 9 ) 

Out of equilibrium the fluctuation-dissipation relation 
Eq. Q does not hold. However within linear response in 
the applied bias Ayt (we consider two leads and without 
loss of generality A/r > 0), we can write an expression 
relating the noise correlator and the dissipative matrix 
7 that generalizes Eq. Q in the limit of low frequencies 
(as compared to the characteristic energy scales of the 
quantum environment, this statement will be made more 
precise in the following sections). We state here the result 
which will be proven later in the text 



( 10 ) 


where we defined uj± = wi A/i, and all remaining quanti¬ 
ties are evaluated at zero frequency and to linear order in 
the applied bias. In general the noise correlator depends 
on both temperature and bias, which at zero frequency 
we denote by Djt’a/x]- It is easy to see that for Ayi = 0 
Eq. (10) reduces to the equilibrium identity ([^, and in 
particular = Dd. For zero temperature but finite 

bias accordingly we obtain 


D(a;) = I + 1^1 ’ 1^1 « 

1 , |w| » Am 

( 11 ) 

The explicit expression for D[o,A/i] is given later in the 
text in Eq. ( [4^ .1^ One should note the similarities and 
subtle differences between the finite temperature expres¬ 
sions, Eqs. 0 , 0 , and finite bias Eq. 0 in the respec¬ 
tive limits of high and low bias and temperatures when 
identifying Ap, ^ 2 T . All these expressions are eval¬ 
uated to first order in ui. From this we note that the 
first correction to Eq. 0 is of order 0{lo^). That the 
linear term in uj vanishes in equilibrium is easy to see by 
looking at the corresponding finite bias result and iden¬ 
tifying D[o,A/.t]AM ^ D[r_o]/ 2 T = 7 ®'! in the first line of 
Eq. ( [ll| ). 

In the next section we show that the Loschmidt echo, 
within a perturbative expansion in the change of the scat¬ 
tering potential, is directly related to the noise correlator 
of Eq. 0. 
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III. LOSCHMIDT ECHO: GENERAL RESULTS 
FOR SMALL DISPLACEMENTS 

The Loschmidt echo in quantum systems is given by 
the (squared) overlap of eigenstates of an initial sys¬ 
tem which evolved in time with two different many-body 
Hamiltonians. Alternatively, it can be seen as measure 
of how close to a given initial state a system comes back 
to, when the evolution on the time reversed path is de¬ 
termined by a different Hamiltonian from the forward 
evolution one. This can be generalized for initial states 
which are a quantum statistical mixture.!^ Denoting the 
Loschmidt echo by the function £(t), this is given by 
£(r) = |A(t)|2, with 

A{t) = ( ) , (12) 

where (...) is the quantum statistical average character¬ 
ized by the initial Hamiltonian "Hi {h = 1) and T-Lt de¬ 
notes the perturbed Hamiltonian. The overlap A{t) is 
denominated fidelity amplitude. 

The relation between the Loschmidt echo and the or¬ 
thogonality catastrophe is seen by treating Anderson’s 
orthogonality as a dynamical process.^ In the problem 
of X-ray absorption spectrum of a metal, the creation of 
a deep hole produces a “shake up” of the Fermi sea that 
causes a suppression of Mahan’s powerlaw divergence at 
treshold frequency (known as the Fermi edge or X-ray 
singularity),^ with an exponent that can be identified di¬ 
rectly with Anderson’s orthogonality exponent .1^^ This is 
captured by the hole propagator which can be calculated 
by evaluating the overlap of the fermionic ground state 
evolved with a Hamiltonian including the core hole, with 
that of the ground state evolved without the hole. This is 
therefore nothing else than the Loschmidt echo where the 
two Hamiltonians Ffi, Hf correspond to considering the 
system with or without the potential of the core hole. 
Beyond the original problem of the Fermi edge singu¬ 
larity in metals, different problems in which some local, 
time varying degree of freedom interacts with a fermionic 
environment, can be treated with the same methodol¬ 
ogy and hence some incarnation of the fidelity ampli¬ 
tude and Loschmidt echo appears naturally. Examples 
include the absorption spectrum of Luttinger liquidpSSll 
and beyond,!^ the single channel Kondo problercPH or 
time-dependent impurites in cold atom systems.l^ 

In order to make the connection with the Langevin 
equation discussed above, in this section we obtain an 
expression for the Loschmidt echo via a perturbative ex¬ 
pansion for small changes in the potential Vx, by consid¬ 
ering Hi =^.0 + Vx and "Hf = 'Ho + Vx-i- 5 x, cf. Fig. 
with (5Hx = Vx+ 5 x — Vx small with respect to Hi- An 
important factor to determine the time dependence of the 
Loschmidt echo, is how rapidly the change in the coupling 
potential occurs. This rapidity is determined by what is 
called the “quench protocol”. Here we consider this is 


symbol 

Sudden (P = S') 

Adiabatic (P = A) 

Xp 

1 

2 

ap 

2 

1 

Up 

1 

1/3 

dp 

2 

1/2 


TABLE I. Protocol-dependent constants used thronghont the 
text. Note that while aA is universal and valid for any adia¬ 
batic quench, the rest of the “adiabatic” constants are valid 
exclusively for the linear quench — they are actually depen¬ 
dent on the adiabatic quench protocol. 


given externally by an arbitrary function g(t) such that 


H(t) = Ho + Vx + g{t) S'H:x. 


(13) 


Initially, ( 7 ( 0 ) = 0 so that we obtain Hi- We impose the 
quench is completed at some time r by setting g(r) = 1. 
We consider an open quantum system in which the elec¬ 
trons spend on average some finite time td in the scat¬ 
tering region, is refered to as the dwell time and we 
consider it to be the smallest time scale in the system, in 
the spirit of the Born-Oppenheimer approximation. This 
defines the time scale for the quench. In particular in 
this work we will study two complementary protocols: 
sudden and adiabatic quenches. In the sudden quench, 
the scattering potential is changed suddenly, which is re¬ 
alized by a step-function shape such that g(t) = 1 for 
t > 0. For the adiabatic quench instead, the potential is 
ramped up slowly, where slow refers to the dwell time. 
For the adiabatic quench, we choose a linear ramping pro¬ 
tocol g{t) = ijr. We show in Appendix [A| that, in the 
limit T ^ Tjo, our results in equilibrium are independent 
of this choice, while for an imposed bias our results are 
characterized by coefficients that depend on the specifics 
of the adiabatic protocol. In what follows we will treat 
in parallel both abovementioned protocols, we list in Ta¬ 
ble ITH] the protocol-dependent parameters as used in the 


text. The Hamiltonian defined in Eq. (13) is time de¬ 


pendent through g{t). To treat this time dependence we 
write the fidelity amplitude in terms of evolution opera¬ 
tors 


A[T) = {ul{r,Q)U{r,Q)). 


(14) 


The operator Uq is the time-evolution operator of the 
(constant) initial Hamiltonian Ho + Vx and U that of 
the (time-dependent) Hamiltonian H(t). In the case of a 
sudden quench, we recover the usual expression Eq. (121. 


We introduce now the interaction picture with respect to 
the initial Hamiltonian. This allows us to write 


^(r) = (T exp - z J dt g{t) ) (15) 

0 


with time-ordering operator T and (5Hx(f) = 

gi«ii^Hxe“*^‘*. The expression given in Eq. ( [T5| ) is the 
starting point for the perturbative expansion given be- 
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low. 

We now perform a pertubative expansion of the fidelity 
amplitude in the displacement SX up to the first non-zero 
terms both in imaginary and real parts. We assume the 
scattering potential to be well behaved, such that small 
changes in X correspond to small changes in V.^ The 
corresponding change in the Hamiltonian is 

S?ix = ^ daVx^Xa = ^ daJixSXa ■ (16) 

a CK 

We therefore obtain 

r 

In^(T) = -i J dtg{t) {Sn:^{t)) 

0 

T r 

(17) 

where Dap{t,t') is the noise correlator given in Eq. ([^ 
and we have used Eqs. ^ and Q. As anticipated in 
the previous section, due to the sum over the indices a 
and /3, only the symmetric component of the noise cor¬ 
relator is relevant. The second term in Eq. p7| ) is a real 
quantity, while the first order term is purely imaginary 
and hence contributes as an overall phase. This phase 
is directly related to the infinitesimal work made by the 
Born-Oppenheimer force, which is consistent with the 
shift in the dynamical phase of the system’s eigenstates 
which is acquired due to the change in the potential SHx- 
This can be seen by expressing the quantum statistical 
average {dHx) in terms of scattering states.l^ For the 
adiabatic {P = A) and sudden (P = S) quenches we 
obtain 

Ap{t) = e (18) 

with Aa = 2 and A 5 = 1 , where the reduction by a factor 
of 2 of the adiabatic fidelity phase with respect to the 
sudden quench results from integrating over the linear 
ramp-up of the potential. 

The Loschmidt echo, in turn, is given solely in terms 
of the integrated two-time noise correlation function 

r r 

ln/:(r) = - J dtj dt'git) g{t') 5X^ ■ D(t, t') ■ 5X , 

0 0 

(19) 


where the function g{t) enters as a weight factor. We 
note that assuming Gaussian white noise, where D(t,t') 
is delta-correlated in time, it immediately follows from 
Eq. (19) that the Loschmidt echo decays exponentially 
with a strength proportional to • D(X) • i5X in the 
large time limit. For now we keep the results general 
and discuss the regime of applicability for the white noise 


limit later. Expressing the fluctuating force correlator by 
its Fourier transform we readily observe that the decay 
of the Loschmidt echo is determined by the symmetric 
noise correlator D(a;),l^ 


, , f duj 

k£(r) = - /- 


dtg{t) i 


(5X^ • D(a;) • (5X. 


( 20 ) 


The time integral above can be performed once the 
quench dynamics g(t) is specified 


00 

\nCp{T) = -J ‘^Bp{u;,t)SX'^ ■-D{uj)-SX, 


( 21 ) 


where the function Bpiui^r) is protocol-dependent 

D £ ^ o 1 - cos(a;T) 

■Bs(w, r) =2-5- 

^ ( 22 ) 
, , 2 [1 — cos(a;T) — wr sin(a;r)] -I- 

Pa(w, t) =■- _2, .4 ---- ■ 


The dwell time Tp) is a characteristic time scale for the 
scattering of the fast (electronic) degrees of freedom and 
provides a high energy cutoff l/rp for the energy inte¬ 
grals in Eq. (211. At the same time, the function Bp{uj, t) 
selects frequencies w < 1/r. This allows us, in the limit 
of large t :$> Tp, to neglect the dynamics of the fast de¬ 
grees of freedom and evaluate the fluctuating force cor¬ 
relator D(a;) in Eq. (211 in the limit of small frequencies, 
uj ^ 1 /t^O — we will use this fact in the next sections. 

The expressions obtained in this section are, within 
the limit of validity of the perturbative approach, quite 
general. In particular, they hold for out-of-equilibrium 
situations. To investigate these results we start in the 
next section with the equilibrium case, for which we can 
straightforwardly apply the fluctuation-dissipation theo¬ 
rem as given in Eq. ( [7|). The out-of-equilibrium regime 
is treated later in Sec. IVI 


IV. EQUILIBRIUM 


We consider here the equilibrium case for which all 
leads are kept at a same chemical potential denoted by g. 
In equilibrium the Anderson orthogonality exponent has 
been shown in Ref. HD to be proportional to the friction 
coefficient of the noninteracting fermionic environment 
for finite systems. This corresponds to the t —>■ 00 limit 
of the Loschmidt echo.^^ We generalize this result here 
to the case of an open system with a continuous energy 
spectrum by calculating the decay of the Loschmidt echo 
for finite times r. 


We continue with evaluating Eq. (21) at T = 0 and 
we discuss the effect of finite temperature later. We can 
therefore use the fluctuation-dissipation theorem as given 
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in Eq. We conclude that, to leading order in 

ln/:p(r) = -i J ujdujBp{u},T)SX^ (23) 
0 

where 7 “'^ is the equilibrium friction coflicient evaluated 
at zero frequency. We see therefore that, as expected, 
in equilibrium the Loschmidt echo is closely related to 
dissipation. We obtain 


In £p(t) = — 


O' p 


In ( — j + 7e 


,5X^ • 7®'^(X) • 5X 


(24) 


where ap is protocol dependent, with as = 2 , = 1 

for the sudden and adiabatic quenches respectively, and 
7 e = 0.5772 is the Euler-Mascheroni constant. We show 
in App. 1^ that the value aA = 1 is independent of the 
assumption of linearity for the adiabatic quench protocol. 
Therefore, 






(25) 


so that the decay of the Loschmidt echo in equilibrium, 
both in the sudden and adiabatic quench scenarios, is 
a powerlaw controlled by the friction coefficient of the 
fermionic system. The powerlaw decay of the Loschm idt 
echo is consistent with known literature result^^^^ and 
reflects the Anderson orthogonality catastrophe.l^^ Eor a 
finite system the r —>■ 00 limit of the Loschmidt echo 
can be obtained, up to prefactors, by replacing r/rp by 
the number of particles of the system and the powerlaw 
takes the usual Anderson’s form. This is justihed since 
the ratio r/rp can be taken as an estimate of how many 
particles have been scattered up to time r — for r —>■ 
00 , all particles in the system have participated in the 
scattering. 

Note that the powerlaw decay of the Loschmidt echo is 
inconsistent with a delta-correlated noise — recall that 
white noise implies an exponential decay of the Loschmidt 
echo. The powerlaw decay signals the breakdown of the 
Markovian, semi-classical Langevin equation Q in equi¬ 
librium and at zero temperature, for which case the clas¬ 
sical noise correlator is zero. In other words, the system 
loses its memory as a powerlaw in time instead of expo¬ 
nentially, which renders the Markovian approximation in¬ 
applicable. An exponential decay of the Loschmidt echo 
is recovered either by imposing finite temperature or a 
finite bias voltage. In the following we comment on the 
finite temperature case, and we reserve the next section 
for out-of equilibrium effects. 

Eor temperatures such that T 3> l/r, we can use the 
classical version of the fluctuation-dissipation theorem 
Eq. (|^ in Eq. (19) to obtain 


and therefore we recover, for r ^ rp, an exponential 
decay governed by the thermal noise 


Cp{t) = e 


_ ^-tPp 5Xt-D 


[T,0] 


■<5X 


(27) 


where j3p is protocol dependent, with /3s = 1, /3^ = 1/3 
for the sudden and adiabatic (linear) quench respectively. 
Note that in this case the coefficient Pa depends on the 
nature of the adiabatic protocol. 

As a last remark of this section, we observe that in 
equilibrium and zero temperature the adiabatic and sud¬ 
den Loschmidt echo are related by a simple exponent .1^ 
Eor zero temperature, from Eq. (24) we obtain 


Cs{T)=CA{Tf 


(28) 


As pointed out before, ua = 1 is independent of the adi¬ 
abatic protocol, and therefore Eq. (28) holds generally. 


The relation given by Eq. (28) has been recently pointed 
out for particular examples in Refs. |43] and |44] for infi¬ 
nite T in finite systems, and argued to be valid in more 
general situations.!^ Eor hnite temperatures we obtain 
instead from Eq. 


Cs{T)=CA{Tf/^\ (29) 

which is valid to leading order in Tpijr. 

V. OUT-OF-EQUILIBRIUM 

In this section we take a step further and allow for the 
presence of an applied bias voltage, represented by dif¬ 
ferent chemical potentials in the leads. Eor clarity we 
consider only two leads which are kept at a chemical 
potential difference A/x > 0. We obtain the decay of 
the Loschmidt echo in the limit of linear response, for 
which the applied bias is small as defined by the condi¬ 
tion AfiTp ^ 1. By evaluating Eq. (21) with the out- 


of-equilibrium noise correlator given in Eq. (10), we can 


obtain closed expressions for the time dependence of the 
Loschmidt echo in terms of the macroscopic coefficients 
appearing in the Langevin Eq. ([^. These expressions 
are valid for all times longer than the dwell time, as de¬ 
tailed below. The derivation is lenghty and therefore we 
summarize here the main results and give a sketch of the 
calculation in the next subsections, while the details can 
be found in the corresponding appendices as listed. 

It is instructive to consider the long- and short-time 
dynamics of the Loschmidt echo as compared with the 
timescale determined by the inverse of the imposed bias, 
since in these limits the expressions simplify considerably. 
We state here the results for zero temperature. The case 
of short-time dynamics (but still large times compared 
with the dwell time) is given by the condition A/xr 1. 
Eor these short times the system is being probed at high 
energies and it is not sensitive to the applied bias. We 


In Cp (t) =-2 Ppt T5X^ ■ ■ 6X , 


(26) 
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therefore recover the equilibrium result 
£p(t) oc ( — 

\td 


(30) 


Here 7 is the full dissipation matrix evaluated to first 
order in the bias. Ignoring this first order correction due 
to the bias we recover exactly the equilibrium result ob¬ 
tained previously in Eq. (24). 


In the opposite limit of long-time dynamics, A fir ^ I, 
the system is more sensitive to the non-equilibrium im¬ 
posed by the bias which results in a different qualitative 
behavior. The major effect due to bringing the system 
out of equilibrium is an exponential suppression of the 
Loschmidt echo in the long-time dynamics, compared 
with the equilibrium power-law decay in Eq. (24). We 
obtain 


£p(t) oce (A/rrp))" 




t ^[O.Ari 


An 


^ -I^^|.5X 

td) 


•iX 


(31) 


where the dissipation matrix and noise correlators are 
evaluated to first order in the bias. The exponential sup¬ 
pression of the Loschmidt echo is dictated by the shot- 
noise fluctuations in the system given by the noise cor¬ 
relator D[o,Ap] to first order in the bias.l^ Note that 
this exponential decay is completely analogous to that 
for equilibrium and finite temperatures in the long-time 
limit (Tr 3 > 1) given in Eq. (27), which is dominated by 
the termal noise. 


The exponential decay in Eq. (31) comes on top of a 


powerlaw behavior, with an exponent that is also mod¬ 
ified from equilibrium showing a competition between 
fluctuations and dissipation, and with the possibility of 
a change of sign in the exponent. This powerlaw correc¬ 
tion is absent in the complementary case of equilibrium 
and finite temperatures given in Eq. (27), due to the van¬ 


ishing of the linear order in frequency correction for the 
thermal noise correlator in Eq. ([^. Note that the out- 
of-equilibrium powerlaw crosses over to the equilibrium 
one at a time r « 1/A/r as expected. 

We proceed now with the derivation of these results. 


A. Out-of-equilibrium noise correlator 


We see from Eq. (19) that the force-force noise correla¬ 


tor is crucial to determine the behavior of the Loschmidt 
echo. In this section we derive the expression given in 
Eq. (10) for the noise correlator within a scattering ap¬ 


proach which highlights the connection to the current- 
induced forces in the Langevin Eq. 0- An equivalent 
derivation in terms of Keldysh Green’s functions is given, 
for completeness, in App. Alternatively, the correla¬ 
tor can be calculated within a Feynman Vernon influence 


functional approach.l^ 

We now proceed with evaluating D(t, t') as given in 
Eq. 0 in terms of single-particle scattering states. For 
this we introduce the notation 


dc.Vi^ie,e') = {^pf+{e)\d^V^\^p^+is')) (32) 


for the matrix elements of the rmresentation of VVx in 
the scattering basis — cf. App. ^ Here |'0^+(e)) is the 
single-particle retarded scattering state with combined 
channel-lead index n and energy e. For notational con¬ 
venience in what follows, we further define the function 

= {d^V^-[e,e')dpV^\e',e)]^ . (33) 

Using Eq. Q to evaluate the quantum statistical ex¬ 
pectation valued^ appearing in the noise correlator 0 
we obtain 


kn 

X (£,£')• 


(34) 


After Fourier transforming we obtain a general expres¬ 
sion for the force-force noise correlator as a function of 
frequencies 

i>H=/gi:A(-i) 

kn 

/ a; uj\ 

X Kfe„ ' 

We observe that the function Kfc„(£, e') contains overlaps 
of scattering states which are are associated with scatter¬ 
ing events including an energy transfer lo = e — e'. We 
expect these overlaps to vary within energies up to the 
inverse dwell time I/tjj. We can expect our description 
of the Loschmidt echo in terms of scattering states to be 
valid in the limit t tjo — a description on microscopic 
time scales smaller than td is beyond an adiabatic scat¬ 
tering formulation. Hence in the following we restrict 
our calculations to this limit, and evaluate (£,£') to 
first order in uj ^ 1/td. We will see later in the text 
that this is enough to capture the leading behavior of 
the Loschmidt echo as a function of time. It is also con¬ 
venient to define the function 



Kfe„(£) = Kfc„(£,£). (36) 

We note that the function Kfe„(£) is closely related to the 
dissipation matrix in equilibrium. At zero temperature 
and taking all leads to be at an equal chemical potential 
jjL, as we show in Appendix it takes the simple form 

^^K,„(/r)=7‘^^ (37) 

kn 

in agreement with the result found previously by use of 















the fluctuation-dissipation theorem. 


The second contribution to Eq. (35) is given by the 


product of Fermi functions /fc(e)[l — /„(e')]. Due to this 
product, the average energy e = {e + e')/2 is limited to 
a region of size A^kn = — IJ-n around the respective 

average chemical potential = l/2(/ifc To make 

analytical progress, we limit our results to the linear re¬ 
sponse regime Afikn td ^ 1 , which allows a perturbative 
treatment of the function Kfc„(e,£') for small deviations 
of e around 

Given these considerations, we calculate Eq. ( [3^ to 
leading order in tjj/t, in the linear response regime (lin¬ 
ear order in A/iru). Expanding Kfe„(£,e') to first order 
in Lu and e — we obtain 

Kfen(e,e ) = 

(38) 


where we have introduced 


dl^^'Kknix) = -{de± de') Kfc„(£,e') 


(39) 


derivatives of Kfe„. We note in passing the following use¬ 
ful properties: Kknie) = K„fe(£), d^Kknie) = 9|K„fc(£) 
and d^Kknjs) = Substituting Eq. (38) into 

Eq. (351, the energy integral can be performed to obtain 


dm= 5 -i: 


W -I- A^lkr. 


which describes the symmetric and antisymmetric energy 


27 T g(u;+A/ifcri)/T’ — ^ 

kn 

X [KU-p.J-cud^KU-pkJ] . 

(40) 

To be consistent with the linear response approximation, 
the functions Kfc„ and have to be evaluated to 

first order in A^td. We proceed with this expansion 
below for the case of two leads, for which the expressions 
are more transparent. The indices k, n describe hereafter 
(the two) lead indices only, where we implicitly assume a 
summation over the channel index. A generalization to 
an arbitrary number of leads is straightforward. Without 
loss of generality we write = /r — A/i/2, /i^ = ^ + 
A/i/2 with A/i > 0. Explicitly, = Kii(/i) -|- 

Afid^KLLifJ.) and Knnifin) = Ki^i^(/i) - A/i9|Kflfl(/i). 
Hence, within linear response we obtain (with u)± = w ± 
A/i) 


D(cc) 


^ w coth [KLL(/i) -t- Kijii(/i) -1- A/i {^l'KLL{^J) - dl'Kimip))] 

coth (^) + w- coth KLij(/i) 

-w+ coth + uj- coth d^KLR(fi). 


1 

dir L 

OJ 

dir . 


(41) 

(42) 

(43) 


The coefficients appearing in the expansion can be in¬ 
terpreted in terms of the different dissipative contribu¬ 
tions with the help of the relations found in App. 
The connection to the friction tensor is found by expand¬ 
ing the dissipation tensor 7 to first order in A/itd a^^ 
7 = 70 ^^ -I- 7 ®'^ -I- 7 ”®'^ -I- ..., where the subscript 0 (1) 
denotes the zeroth (first) order in the expansion respec¬ 
tively of the equilibrium (eq) and non-equilibrium (neq) 
contributions to the friction tensor. In App. 0 we show 
the following identities 


ir = ^[dt^RL{k^)-diKM] 

D[o,A/.j] = [KLfl(/r) -f K_RL(Ai)] , 


which together with Eq. (371 imply 


KLL{^k) + ^RR{^k) = 4^ ( 7^ - 


(44) 

(45) 

(46) 

(47) 


r 


Plugging in these identities in Eq. (41) we obtain the 
anticipated result stated in Eq. (10). 


In order to be sensitive to non-equilibrium effects, we 
need to impose temperatures smaller than the bias. We 
therefore take the zero temperature limit of the noise 


correlator given in Eq. (10). For |a;| < A/i we then obtain 

D[ 0 .AAt] 


D(w) = D[o,am] 


w 


7 


eq 


A/i 


a ;2 7“q 

T A/i 


while for |w| > A/i we have 

D(w) = |w|7, 

where we have used that lim 3 ;_,,j-oo cothx = ± 1 . 

B. Short- and long-time dynamics 


(48) 


(49) 


Inserting Eqs. (48) - (49) into Eq. (21) we can express 


the decay of the Loschmidt echo in terms of the meso- 
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scopic coefficients that control the Langevin dynamics of 
a heavy particle embedded in the fermionic environment, 
cf. Eq. 0 . We obtain 


ln/:p(T) = - ^ 
TT 


7 . + ln|^) 


Qp 

TT 


7 e + In (A/xr) 




(5X^ • 7 • 5X 

D[o,Am] 


Afj. 


+7ri 


(5p 


- /3p T 5Xt . D[o,Ap] . <5X - ^ (5Xt . 7f• 5X 


+ — fap — 2/3p cos(A/xt)') (5X^ • —. 

TT V / A// 


^X. 


(50) 


with (5s = 2, (5 a = 1/2. Eq. (50) gives the behavior of the 
Loschmidt echo at arbitrary times larger than the dwell 
time, up to quadratic order in (5X. In the following we 
further investigate different timescale regimes. 

The short-time dynamics is given by the limit Afj, r <C 
1 — “short times” here should be considered as short 
with respect to the inverse bias time scale but long com¬ 
pared to the dwell time rp. In this regime we conclude 


ln£p(T) = - — 

TT 


7e ' 




^Xl • 7 • (5X, (51) 


which yields Eq. (30). Note that the equilibrium friction 


term given by 71 /^ constitutes the dominant contribution 
for the decay, and the full friction matrix 7 is restricted 
to positive values for small A/xrp ^ 1 , which ensures 
that Cpq(t) < 1 for all times r. 

We turn now to evaluating Eq. ( |50| ) in the long-time 
limit A/xr ^ 1. Writing ln(T/rp) = ln(A/xr/(A^rp)) 
and observing that A/xrp In(A^rp) goes to zero for 
A/xrp ^ 1 we obtain 


ln£p(r) = -/3pr5Xt.D[o.AH.^] 
— — 7 e -f In (A/xr) (5X1 . 

-I- — In (A^rp) (5X1 . . jx 

TT 


Dr 


7 ^Xl • 7 f‘1 • (5X 

TT 

D 


— ^ap — 2f3p cos(A/iT)^ 5X1 


from which we obtain Eq. (311 by keeping the dominant 
terms in large r.l^ 

The exponential suppression of the Loschmidt echo 
is dictated by the shot-noise fluctuations in the system 
(note that D[q A^i] is positive definite) and it is consis¬ 
tent with a leading behavior of Gaussian white noise for 
the fluctuating force correlator D(t,t') in Eq. (0. Fur¬ 
thermore, the exponent of the powerlaw shows a com¬ 
petition between fluctuations and dissipation, which is 
a clear signature of the departure from equilibrium — 


as remarked before, the corresponding powerlaw is ab¬ 
sent in the equilibrium, finite temperature case. Since 
7 ^^ = D[p_q]/( 2 T) for zero bias, the sign of the expo¬ 
nent —5X1 . _ D[o,a/j]/(A/x)] • 5X depends on the 

assymetry between shot and Nyquist noise in the lin¬ 
ear response regime. In fact, this exponent can be pos¬ 
itive for finite bias, which leads to an enhancement of 
the powerlaw instead of the usual decay, which has been 
dubbed as “anti-orthogonality” We note however that 
—5X1 _ D[o a ui/(A/x)] • 5X is always negative (see 

Eqs. ( [47| ) and ( |C.4[ )). The sign change of the exponent 
happens when the leading order Kpp(/i) « —Kpp(/x) 
cancels out, such that the linear order correction 7 ®'^ in 
the bias becomes dominant, which can then lead to a 
change of the sign of the exponent. 

Since the out-of-equilibrium short-time dynamics is 
essentially the equilibrium one, the identity £s(t) = 
CAirY is still fullfilled in this limit, as can be seen di¬ 
rectly from Eq. (51). On the other hand, in the large¬ 
time regime, we conclude from Eq. (52) that this iden¬ 
tity is violated due to the factor j3p in the exponential. 
This difference, attributed to the structure of g{t), can 
already be obtained by looking at Eq. (19). Assuming 
white noise, we immediately deduce from Eq. (19) an 
exponential decay of the Loschmidt echo with an expo¬ 
nent —/3p5Xl • D[o_A/i] • 5 Xt. The powerlaw decay in 
Eq. (52) constitutes minor correction terms to the white 
noise assumption, and therefore to leading order in rp/r 
we obtain £s(t) = as in the equilibrium, fi¬ 

nite temperature case given in Eq. (29). 

In Sec.|Vl]we study these results for a specific example 
of a two-level model coupled to one vibrational mode. 


VI. EXAMPLE: TWO-LEVEL MODEL WITH 
ONE VIBRATIONAL MODE 

In this section we analyze the sudden quench 
Loschmidt echo for the example of a system with one 
classical degree of freedom connected to two leads. This 
serves as an toy model to illustrate the above results. 
The “heavy” classical degree of freedom X = X{t) corre¬ 
sponds to a mechanical vibrational mode of the system. 
Accordingly we consider the Hamiltonian 

H = nx + nL + nD + nT (53) 
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where the different terms are specified as 

'Hl= [ ^^ri)cl{e)cn{e) 

dm' 

mm' 

X! {4j(.£Wv'n{£)dm + h.C.) . 


(54) 

(55) 

(56) 

(57) 


Here, the operator cl^{e) [c,,(e)] creates [annihilates] elec¬ 
tronic states \(j)jj{s)), which are approaching the scat¬ 
tering region from lead rj = L, R with chemical poten¬ 
tial > Hr- Hx describes the evolution of the pa¬ 
rameter X with potential U{X), mass M and frequency 
ujQ. Hr models the two-level system (quantum dot) 
with states jm), created (annihilated) by the operators 
din (dm)- Ht represents tunneling between the leads 
and the system with tunneling amplitudes Wrjmis) = 
{(l)n{s)\W\m)/y/^- The coupling of the mechanical de¬ 
gree of freedom and the electrons in the dot is described 
by the matrix ho{X). 

We consider a two-level system with degenerate energy 
levels So- The single oscillator mode X is assumed to 
couple to the difference in the energy level occupation 
with a strength given by A. Hence we write 


ho{X) — 


e + t 
t E- 


(58) 


with interdot tunneling amplitude t and e± = eo i AX. 
Tunneling from the left (right) lead to the two-level sys¬ 
tem and back is described by the amplitudes T^ (T/j) 
which for simplicity we take as Ti, = T^j = r/2. In the 
wide-band approximation these amplitudes are assumed 
to be energy independent. With these definitions, the 
coupling matrix W^jm reads 

_ ( v/tTM 0 \ , 

*'-( 0 vWfS))’ ' ’ 

and the frozen retarded Green’s function takes the form 


Gxi^) — 


Ax(£) 


e — E- + iT/2 t 

t 6 — (S_i_-i-ir/2 


(60) 


with Ax(e) = (e —e_ -l-ir/2) {e — e+ +iT/2) This 
model was studied in Ref. [5^ in the context of current- 
induced forces. The frozen scattering matrix S and its 
first order non-adiabatic correction A are given by 


Sx[e) = 1 - 


iT 

Tx(e) 


1 1 
1 1 


(61) 



FIG. 2. Sudden quench Loschmidt echo as given in 
Eq. ( |50[) and comparison to the short- and long-time behavior 
Eqs. ( |30| ,( 311 (with the corresponding proportionality coefli- 
cients), for the example of a two-level system coupled to one 
vibrational mode with T = t = 0.48, td = l/F = 2.08, eo = 0, 
Ar = 0.04, SX = 0.15, X = 0.4, T = 0. All energies (and 
inverse times) are in units of A^/ (MuJq) and distances in units 
of A/(Mcn§). 


XsM 


Xs(t) 




(a) 


(b) 


FIG. 3. Linear response solution of the sndden Loschmidt 
echo — Eq. ( |50[ ), vs. general solution [i.e. Eq. ( |35| l inserted 
into Eq. (211] at different bias voltages, (a) Ar — 0.04, (b) 
Ar — 0.4; others parameters as in Fig. 


^x(£) = 


An 


0 1 
-1 0 


(62) 


where Lx{e) = e — e+ + iT . With the expression of 
the S-matrix and the A-matrix, we can determine all 
the mesoscopic coefficients appearing in the Langevin 
Eq. ([^, which determine the behavior of the Loschmidt 
echo. This is depicted for arbitrary times t > tr in 
Fig. [^according to Eq. (501 , and compared to the small 
and large time regimes expressions in Eqs. (30),([M|). The 
dwell time enters in this example as a timescale which is 
of the order of the inverse tunneling amplitudes. 

The considered model allows us to analyze the 
Loschmidt echo also outside of the linear response regime, 
by directly evaluating the colored noise-noise correlator 
as given in Eq. (35 1 . The matrix 9a^6c(£, ^0 [see Eq. (32 1 ] 
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E(X) 



FIG. 4. (a) The out-of-equilibrium powerlaw coefficient 

E(X) = —^( 70 '* -I- qj'* — for A/i = 0.04 as a function 
of A; other parameters as in Fig.j^ (b) Zoom-in of (a) which 
shows the change of sign in the exponent. 


is given by 


doyl^{e,e') = 27 r[W- G§(eydxho(X)G§(e') ■ . 

(63) 


Eq. (35) is neither restricted to A/ttd ^ 1 nor to the 


regime T^jT ^ 1 and hence is valid for arbitrary A/r 
and all r. Thus we can study the Loschmidt echo for 
increasing bias voltages A/x. A comparison of the general 
solution, that is substituting Eq. (35) into Eq. (21), and 
the linear response solution Eq. (50), is depicted in Fig. 
The figure shows that the linear response solution agrees 
very well with the general solution for small A/x tjj . 

To close this section, we show that the powerlaw expo¬ 
nent — ^( 7 ®'^ — D[o,A/ij/A^) can present changes in sign. 
This is shown in Fig. |^for a specific finite bias voltage, 
where the exponent becomes positive for a small range 
of displacements. 


VII. SUMMARY 


sipation coefficient of the heavy particle, in agreement 
with the results found in Ref. [21] in which the exponent 
of the Anderson orthogonality was related to the dissipa¬ 
tion of a heavy particle moving in a finite-sized quantum 
environment. For finite temperatures, in the limit of long 
times we recovered an exponential decay which reflects 
the classical version of the fluctuation-dissipation theo¬ 
rem. 

When a small bias is imposed, we showed that in the 
long-time dynamics (as compared to the energy scale 
given by the bias), the Loschmidt echo is dictated by 
an exponential decay with a strength given by the shot- 
noise fluctuations. The exponential decay is consistent 
with a white-noise spectrum to leading order. As a cor¬ 
rection term to white noise, the Loschmidt echo shows 
an algebraic behavior with a powerlaw exponent given 
by a competition between fluctuations and dissipation. 
This competition can give rise to changes in the sign of 
the powerlaw exponent. The powerlaw correction to the 
exponential decay is characteristic of non-equilibrium ef¬ 
fects and is absent in equilibrium at finite temperatures. 
In the case of short-time dynamics, the system is mostly 
insensitive to the applied bias and the Loschmidt echo 
still presents a powerlaw decay, controlled by the full 
non-equilibrium dissipation coefficient. 

The results summarized above are independent of the 
quench scenario. In particular we studied the comple¬ 
mentary cases of sudden and slow quenches, and showed 
that in equilibrium the sudden quench Loschmidt echo 
is the square of the Loschmidt echo for the slow quench, 
independent of the functional form of the slow quench. 
This relation generalizes the relation found for finite 
quantum systems at infinite times in Ref. |43| and 1441 
where a Luttinger liquid model which is subject to a lin¬ 
ear slow quench was studied. We find that this corre¬ 
spondence breaks down out of equilibrium or for finite 
temperatures. To leading order however we can still es¬ 
tablish a simple relation involving a non-universal expo¬ 
nent. 


In this paper we calculated perturbatively the decay of 
the Loschmidt echo for an open non-interacting fermionic 
system in the presence of an external bias, which is sub¬ 
ject to a scattering potential quench. We expressed our 
results in terms of the mesoscopic quantities describing 
the complementary problem of a heavy particle moving in 
a quantum environment, and showed that the Loschmidt 
echo decay is controlled by the noise correlator of the 
heavy particle. 

This result allowed us to study Anderson orthogonality 
effects for the open quantum system. In the particular 
case of equilibrium and zero temperature, the decay of 
the Loschmidt echo is a power law controlled by the dis- 
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Appendix A ADIABATIC PROTOCOLS 


In this appendix we show that Eq. (25) is independent of the adiabatic quenching protocol. In particular we show 
that for any function g{t) which grows as a powerlaw to its maximal value on a timescale much larger than tu, we 
obtain a powerlaw decay of the Loschmidt echo as in Eq. (241 for ap = 1. We are interested in the equilibrium 
behaviour at zero temperature. Equation ( |23[ ) is valid for any quenching protocol with g{0) = 0 and g{T) = 1. 

-T^ with n > 1 integer we read off 


Equation 

Specifying the quench dynamics to g{t) = ( 


l/rn 

^X^7"'1(X)-JX. (A.l) 

0 

We split the integral with respect to u! into a part w < 1/t and w > 1/r. Since the integral where w < 1/t gives only 
a constant (which depends on the value of n), it is irrelevant at large t/tp ^ 1. We are thus left with 


InCAQiT) = -^ J duju! J dt 


1/td 

5X^ •7®‘i(X)-JX. (A.2) 

1/t 

We perform an integration by parts to evaluate the integral in the absolute square values. Eor large r/rp we 
approximate to leading order 


InCAgir) = -^ J dujuj J dt 



Inserting this into Eq. (A.l) we conclude that 


InCAQir) = 


(A.3) 


(A.4) 


up to an irrelevant constant for all g(t) = (y)" with n > 1 integer. Comparing to the large time behaviour of the 
Loschmidt echo for the sudden quench scenario we readily deduce the validity of Eq. (251 independent of the adiabatic 
protocol. 


Appendix B ADIABATIC SCATTERING THEORY 

In this appendix we outline the elements of scattering theory used in the main text. Eor "Hx = 'Hq + Vx assuming 
non-interacting particles, we can write 


Hx 


de f de' 


= /S 


Xf 


kn 


(£)««(£'): 


(B.l) 


with the single-particle Hamiltonian iLx = Hq -j- Vx. The operators a^{e) and annihilate and create, respec¬ 
tively, the retarded single-particle scattering states of the Hamiltonian Hx with energy e and combined 

channel and lead index m, hence it follows 

[Hx]knie,e') = {^f+{e)\Hx\^^+{e')). (B.2) 

The corresponding advanced scattering states are indicated with the superscript (—), that is \tp^~{e)). These scatter¬ 
ing states are solutions of the time-independent Schrodinger equation at every time t [note that we consider X = X(t)], 
and obey the Lippmann-Schwinger equation. The retarded and advanced scattering states are defined by their bound- 
ary conditions: |V’m ^ \^)) has incoming (outgoing) waves only in channel/lead m, and evolved from the free states 

\4>m{£)) at t ^ Too, which fulfill Hq |(()m(e)) = £ \4>m{£)) with energy e. We normalize the scattering states to unit 
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flux. The overlap between retarded and advanced scattering states then defines the S-matrix as 

S^^{s)2nS{e-e') = (^f-(£)|^^+(e')) • 


(B.3) 


From this definition it follows that S'x(e)<S'x(e) = 1, because of the scattering states’ normalization to unit flux. 

For a Hamiltonian which parametrically depends on time via the classical parameter X, Eq. (B.3) gives the frozen 
S-matrix of the system: the solution of the time-independent scattering problem at time t. The exact S-matrix S{e, e') 
is defined through the overlap of retarded and advanced scattering states solution of the full time-dependent problem. 
Assuming a slowly vary ing parameter X, the exact S-matrix can be expressed as an adiabatic expansion in the 
velocity i^^ pM'^^\29 \ 5i \ 52 \ context, “slowly varying” means that the dynamics of X is much slower than the 

electronic time scales. To first order in the adiabatic expansion we can write, in the Wigner representation, S{e,t) = 
A^{e)Xa- The zeroth order is given by the frozen S-matrix. The first non- adiabatic correction to the 
S- matrix is given bjl^ 


(e)|a„Ex |V’™’^(e)) - ^ {e)\daVyi\d,'ilj^+{e)). 


(B.4) 


with da = djdXa and {e)). Note that throughout the manuscript we work in the Heisenberg 

picture, so that there is no explicit time dependence on the states, due to unitarity, the S and A matrix fulfilP^ 


(V’n {£)\dc.V\'4’k (e)) =idaSnk{e), 


i^n (.£)\daV\d^i^f'^ie)) = -Al^{e) + ^ dASnkie) ■ 


(B.5) 

(B. 6 ) 


These relations will be used in the next appendix. 

We can express the reaction forces in the Langevin Eq. ([^ in terms of Eqs. (B.31 and (B.4), which we list here for 
completeness 1241251291 Born-Oppenheimer force is given by 


Fa{X)= [ ^^/„(e)tr{n„4(£)9„5x(e)} 


(B.7) 


where tr{...} denotes a trace over scattering channels, /n(e) = [exp[(£ — /in)/T -|- 1 ]]“^ is the fermionic distribution 
function in lead n with chemical potential /i„, and da = djdXa, H^ is a projector onto channel n. The two 
contributions to the dissipative force as discussed in the main text are in turn given by 


7a^(X) = -^ j ^aE[/n(£)]tr|n„5„S']t(e)a;35x(e)}^ 

n 

7a^'’(X)=^ f ^/„(e)tr|n„ a„5']^(e)A^(£) - A^^(e)9„5'x(e)] } ■ 

^ j m s 

The white-noise fluctuating force correlator is given by 

DapOq = / ^ E T [4(e)5a^x(e)] ^ -Uk- 4(£) a;35x(£)} ■ 

k^TTl 


(B. 8 ) 

(B.9) 

(B.IO) 


In equilibrium it is connected to the friction coefficient via Dap = 2 ksT 7 ®^, where 7 ®'^ is evaluated in equilibrium. 
This agrees with the fluctuation-dissipation theoerem in Eq. 0 in the classical limit w <C T. 


Appendix C IDENTIFICATION OF FRICTION AND NOISE WITHIN SCATTERING THEORY 


Below we derive, within the framework of scattering theory, Eqs. (37) and (44) - (46) of the main text, which give 


the connection of the Loschmidt echo to fluctuations and dissipation. The friction tensor and the noise correlator are 
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given in Eqs. (B. 8 ) - (B.IO). In linear response these quantities read 


(7r)a/3 = \^daSl{n) ■ dfjSt 

{daSl{e) ■ dfjStie)^- {daSl{£) ■ dpStie) 


I7i - 8^ \ 


a. 


i?i?J S — 


( 7 P)a /3 = ^{ (/r) • (m) “ ' d^St 


RR ) s 


(C.l) 
(C.2) 
} • (C.3) 

J S 


twice inserting the complete set of advanced scattering states 1 = / conclude from Eqs. (32) 


We begin by expressing the function K'^^{e), defined in the main text in Eq. (33), in terms of the S-matrix. By 
: the complete set of advano 

and ([33| that 

^aP I 


K\ 


(e) = {{^f^i£)\daVyi\^^+{e)){'iljf+{e)\dpVyi\^f+{e))'^^ 

= {(a<,4(^) • (4(e) • • 

Because of the unitarity of the scattering matrix it readily follows that 

{dcSlin) ■ dpStifi)} ^ = (7^) 


OL0 


(C.4) 


(C.5) 


kn 


at r —?► 0 by referring to Eq. (C.l). We thus have proven Eq. (37). We continue with the derivation of Eq. (44). 
Hereto we note that dlK^^{e) = and d^K'^^{£) = a|if“^(e). The latter property follows immediately 

from the relation (e) = (e) and due to the symmetric summation in the indices a and /3. Hence we get 


A/r 
47r . 


Y. 9e [Kfje) - Kt{e) 

n—L,R 


= ^3. 

Stt 

= h?) 


(daSj^e) ■ dpStie)'^- (daSl{£)dpSt{£) 


RR 


1 Jap 


(C. 6 ) 


by referring to Eq. (C.2) in the last step which proves Eq. (44). In order to derive the relation in Eq. (45) we first 
observe that (e) can be written in terms of the S- and A-matrix. Again, inserting two complete sets of advanced 

scattering states, we find 

d^K^{e) = \{d,-d,,)Kt^{e,e') 


= 2 {(^eV'?^^(£)l^«^x|V'^+(e))(V'^+(£)|5/3hSc|V’?^+(e)) 

+ (e) 15a Vx I V'n ^ (£) ) {e)\dpVyi\d{s)) 

- (V'f^(£)|5aI6c5j|^/>^+(e))(?/'^+(e)|j dpVyi\i;f+(e))] 

L J J s 

^ {5<,4’"'"(£)5e [^r(£)4’”'(£)] dpSl\e)}^ 

ml 

Y{{daSUe) ■ St{e))^^ [sUe)A^{e))^ 


ml 


' nk 


■ Stie))^^ (Slie) ■ . 


(C.7) 


In order to identify the pure nonequilibrium friction tensor in linear response [cf. Eq. (C.3)], we observe that 
5e-^fen(£) = ~^e^nki^)^ where we stress that d^K^^{£) = 0. Restricting to two leads k,n = L,R we thus con- 
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elude 


A/i 

47r - 


n—L,R 


= ^{ {daSjin) ■ ■ daStin) 

- (daSl{fi) ■ A^ifi) - A^^n) ■ daStin)^^^ I 

= (iDc. • 


LL 


(C. 8 ) 


Here we made use of tr{9Q,S'/Hf — A^^daSt} = 0 to realize that the equilibrium contribution in 7 ”^'^ vanishes. We 
conclude with Eq. (45). Finally, we show Eq. (46). For two leads, we write the correlator of the fluctuating force 
apm 


D, 


w = / ^ E /-(") [1 - • 

l. _ r D 


(C.9) 


k,m—L,R 


Making use of fk{s) [^ — fkis)] = —T 9e/fc(e) = T6{e — fik) for small T, we can express the correlator in linear response 
to first order in the applied bias voltage as 


^a/3(X) = 


T 

2t: 


E 

ik,7n—L,R 


Afi 


d,Kll{e)-d,Kl^Je) 


+ 


A/x 
47r - 




(C.IO) 


since — ^ Z^kjs )- With ^ we identify the equilibrium friction tensor in Eq. (C.l) and its non-equilibrium 

correction in Eq. (^C.2|) in the above expression for the correlator. For T —)■ 0 we then obtain Eq. (46). 


Appendix D ALTERNATIVE KELDYSH APPROACH 


This appendix presents an alternative derivation of the results in the main text by using Keldysh Green’s functions 
technique. We consider the generic Hamiltonian given in Eq. (53) and we define the time-dependent Green functions 
of the dot 


g^^,it,t') = -t0{t-t') {{dr^{t),dUt')}) (D.i) 

g:^^,{t,t') = i0{t'-t){{d^{t),dl,{m (D.2) 

G^m'it,t')=-^{dm{t)dl,it')) (D.3) 

G:^n.'it,t')=i{dl{t')d^{t)) (D.4) 

where the curly bracktes denote the anti-commutator operation. We assume stationary states throughout 

this section so that the above defined Green’s functio ns dep end on the difference of the time arguments. With these 
definitions the noise correlator defined in Eq. § read^^mH 

Dapit-t')= ^e*(^-^')(‘-*')tr{A„g>(e)A^g<(e')}^. (D.5) 


with Aa = 9jv„ho(X). The functions G^is) and represent the Fourier transforms of the greater and lesser 

functions. It is sufficient to evaluate the noise correlator in the adiabatic limit as this already guarantees that the 
fluctuation-dissipation theorem is fulfilled. Hereto we introduce the Fourier transform of the adiabatic lesser and 
greater Green functions (e) and (e) with respect to a frozen configuration X and conclude for the Fourier 
transform of the fluctuating force 

tr{A„ G> (e - |) Ap G< [e + |) }^ • 


(D.6) 
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Analogously we introduce the adiabatic retarded and advanced Green functions. Using Langreth rule, we can express 
the greater and the lesser Green function, respectively as 


G<(e) = G^(e)E<(£:) G^{e) 
G>(e) = G^(e)E>(e) G^(e) 


with greater and lesser self-energies 


E<(£) = iJ2fk{£)wHe)Uk{e)W{e) 

k 

E>(£) = ^(1 - /fc(e)) lUt(e) Ukie) W{e) 


(D.7) 

(D.8) 


(D.9) 

(D.IO) 


where = \(j)k{e)){(j)k{£) \ is a projector onto lead space and the coupling matrix W is defined via the Hamiltonian 

in Eq. (57). 

Plugging Eqs. (D.7) and (D.8) into Eq. (D.6) we can write the noise correlator as 

kn 

where we defined the function {£,£') as 

kl^{£,£') = tr |a„ g«(£) w\£) nfc(£) W(£) G^{£)Kp g^{£') w\£') n„(e') W[£') g^{£') 


(D.ll) 


(D.12) 


The expression in Eq. (D.ll) has the same structure as the expression in the main text in Eq. ( |35[ ). Here we proceed 
with the analogous steps to get to the result for the noise correlator in Eqs. (48) and (48) with all its consequences for 
the Loschmidt echo and the fidelity amplitude. In order to show their equivalence we are thus left with identifying 
the coefficients 




47r 

Afi 


kn 


(f) + 




D 

Afi 


OC^ 


47r 

1 

47r 




(D.13) 

(D.14) 

(D.15) 

(D.16) 


for two le ads u nder sy mm etric summation with respect to the indices a and in linear response at zero temperature 
[cf. Eqs. (^), (M) - (®]. We defined k^^{£) = k^^{£,£). 

We begin witntne identification o f the friction tensor. We take the expression of the friction tensor in terms of the 
adiabatic Green functions from Refs.^^HHl 

7a/3 = y ^tr|A„G>{e) A^i9eG<(£)|^. (D-17) 

We immediately conclude for the friction tensor in equilibrium that 

(7o‘')a/3=/^ E tr{A„G«(£)Wt(^)nfc(/r)lU(/r)G^(/x)A^G«(M)Wt(/i)n„(M)W^(M)G^(M)}^ (D.18) 

^ kn^L,R ^ 

since fk{£) (1 — f kk)) = —Tdefk{£) = 0 for T —»■ 0 a s well as f{fi) = ^ and —def{£) = S(£ — /r). A comparison to the 
dehnition in Eq. (D.12) readily results in Eq. (D.13). 

Next we address Eqs. (D.14) and (D.15). Hereto we write the Fermi functions of the left and right lead as /L/_R(e) = 
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f{e) =F ^ def{e). Keeping only terms linear in A/r we conclude after performing an integration by parts 




e—fj, 


(D.19) 




de 

27T 


+ 


f^)) de +K^^(e,/x) 

[def{e)f (e) + K2^^{e) - kf^^s) - kf^{e) - K^^ie) - K^^ie) + kf^{e) + 


The last terms vanish due to the symmetric summation over a and /3 and since K2„ {£) = K^^ie). This yields 
Eqs. ( |D.14[ ) and ( |D.15| ). 

Finally we identify the delta-correlated noise and prove Eq. (D.16) in linear response. We rely on the expression 

^o/ 3 M= y ^tr{A„G>(e)A^G<(e)}, (D.20) 

obtained in the literaturd^mS] terms of the Green functions of the dot. Using Eqs. (D.7) and (D. 8 ) we immediately 
identify 


£> 0 / 3 ( 0 ;)=/^ fk{s) {I - fnk)) kl^iie) 

1 _r ID 


(D.21) 


kn—L,R 


To linear response we find Eq. (D.16) by using k2^{e) = and the above stated relations for the Fermi 

functions. 


We end this appendix by showing the direct equivalence between the function k2^{e,e') defined in Eq. (D.12) and 
K2^{e,e') defined in Eq. (33), that is 


k 


a/3 




(D.22) 


since daHx = daVx- We note that we can relate the Green function of the dot and the scattering states via the 
Lippmann-Schwinger equatiorP^ 


nz3 K‘+(e)) = G«(e) 


(D.23) 


where denotes a projector onto the space of the dot. With the aid of daHx = Hd dahoCX.) and the Lippmann- 
Schwinger equation projected onto the dot’s space in Eq. (D.23) we have 

= {{M£)\Wie) G^ie) G^ie')W^{e') |</)„(e'))(</'n(£')l Wie') G^{e')dpHx G«(e) wHe) |</>fc(e))}^. 

(D.24) 


Due to symmetric summation under exchanging the indices a and (3 and using that n„(e') = |(/)„(e'))((/)„(e')| we 
conclude that 

(g e') = tr{a„i7x G^[e') W\e') n„(e') W{e')G'^{e') dpHx G«(e) W\e) nfc(e) W{s) 

= k2^{e,e'). (D.25) 

This shows explicitly the equivalence to the scattering states approach. 
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